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I review the current status of combing weak gravitational lensing with depth information 
from redshifts as a direct probe of dark matter and dark energy in the Universe. In 
particular I highlight: (1) The first maximum likelihood measurement of the cosmic shear 
power spectrum, with the COMB017 dataset (Brown et al 2003); (2) A new method for 
mapping the 3-D dark matter distribution from weak shear, and its first application to the 
COMB017 dataset (Taylor et al 2003); (3) A new method for measuring the Dark Energy of 
the Universe using purely the geometry of gravitational lensing, based on cross-correlation 
tomography (Jain &: Taylor 2003). I show that this method can constrain the equation of 
state of the universe and its evolution to a few percent accuracy. 

1 Introduction 

Gravitational lensing provides us with the most direct and cleanest of methods for probing the distribution 
of matter in the universe (Mellier 1999, Bartelmann & Schneider 2001). The lensing effect arises due to 
the scattering of light by perturbations in the metric, stretching and contracting bundles of light rays, 
causing the distortion of background galaxy images. Hence gravitational lensing does not depend on 
any assumptions about the state of the intervening matter. These distortions manifest themselves as 
a shear distortion of the source galaxy image (see e.g. Tyson et al 1990; Kaiser & Squires 1993), or 
a change in the surface number density of source galaxies due to magnification (see e.g. Broadhurst, 
Taylor & Peacock 1995; Fort, Mellier & Dantel-Fort 1997; Taylor et al 1998) and can be used to map 
the two dimensional projected matter distribution of cosmological structure. As the matter content of 
the universe is dominated by non-baryonic and non-luminous matter, gravitational lensing is the most 
accurate method for probing the distribution of this dark matter. 

Weak lensing studies have been carried out for a wide range of galaxy clusters, allowing precision 
measurements of cluster masses and mass distributions (see e.g. Tyson et al 1990, Kaiser & Squires 1993, 
Bonnet et al 1994, Squires et al 1996, Hoekstra et al 1998, Luppino & Kaiser 1997, Gray et al 2002). On 
larger scales, the shear due to large-scale structure has been accurately measured by several groups (see 
e.g. van Waerbeke et al 2001, Hoekstra et al 2002, Bacon et al 2002, Jarvis et al 2003, Brown et al 2003). 

Depth information, from galaxy redshifts, has already been used in weak lensing studies to determine 
the median redshifts of the lens and background populations. This can be a limiting factor in the analy- 
sis, leaving uncertainty in the overall mass normalisation. In this review I describe an application using 
accurate shear and redshift information from the COMBO 17 dataset (Wolf et al, 2003) to estimate the 
first maximum likelihood shear power spectrum analysis (Section 2), including all of the uncertainties 
associated with source distances (see Brown et al 2003). While the shear power contains much important 
information on the statistics of the dark matter distribution and cosmological parameters, a lot of infor- 
mation is projected out. However this is not a necessary step. In Section 3 I outline how the full 3-D dark 
matter distribution can be recovered from shear and redshift information, and make the first application 
to the COMB017 data (see Taylor 2001, Bacon & Taylor 2003, Taylor et al 2003). Finally, in Section 4, 
I describe a new method for using shear and redshifts in a purely geometrical test to accurately measure 
the equation of state of the dark energy in the universe and its evolution (See Jain & Taylor 2003). 

It is clear that the combination of shear and redshift information provides a powerful tool for cosmol- 
ogy, allowing us to not only image the dark matter directly, but also see its evolution over cosmic time. 
Hence future lensing surveys should consider being coupled to photometric redshift surveys, opening up 
new dimensions in gravitational lensing studies. 

1.1 Basic weak lensing equations 

The metric of a perturbed Robertson- Walker universe in the conformal Newtonian gauge is 

ds 2 = -(1 + 2$)d< 2 + a 2 (t)(l - 2$)dr t dr\ (1) 



where $ is the Newtonian potential, a is the cosmological scale factor, and we have assumed a spatially 
fiat universe. The Newtonian potential is related to the matter density field by Poisson's equation 



V 2 $ = 4TTG Pm 6a 2 = ^A^^a" 1 ^ (2) 

where S = Sp m /p m is the matter density perturbation, A# = 1/H m 3000 /i _1 Mpc is the Hubble length, 
and £l m is the present-day mass-density parameter. The lensing potential, (f>, for a source in a spatially 
flat universe at a radial distance r is given by (e.g. Bartelmann & Schneider 2001) 



(f>{r,r0) = 2 J o dr ' $ ( r '' r ' 6 ')' 



(3) 



at an angular position on the sky. Hence the lensing potential is just a weighted radial projection of 
the 3-D gravitational potential. 

In lensing the symmetric, traccfree shear matrix, 7^, which describes the distortion of the lensed 
image, and the magnification, p, which describes the change in area, are observables. The shear matrix 
is 

iii = (didi-\&?i&)<l>, (4) 

where di = r(Sij — 9i9j)Vj = r(Vj — 9id r ) is a dimcnsionlcss, transverse differential operator, and 
d 2 = did 1 is the transverse Laplacian. The indices take the values (1,2), assuming a flat sky. The 
lens convergence, k, is defined by 

n= 1 -d% (5) 

An estimate of the lensing potential, 0, can be found from the shear field by the Kaiser-Squires (1993) 
relation; 

^(r,r6) = 2d-%d jltJ (r,rO), (6) 

where d~ 2 = J d 2 9' In \0 — 0'\/2ir is the inverse 2-D Laplacian operator. In practice the shear field is only 
discretely sampled by galaxies so we must smooth the shear field to perform the differentiation. This 
also serves to make the uncertainty on the measured shear field finite, since each source galaxy has an 
unknown intrinsic ellipticity. The scalar convergence field, k, can be estimated from the shear field by 
combining equations (6) and (5), up to a constant of integration. 

In principle there is another odd-party pseudoscalar which can be formed from a shear field; 

(3 = e?d j d m d-' 1 lij , (7) 

where is the 2-D antisymmetric Levi-Civita tensor. This field cannot correspond to the even-parity 
lens convergence field and can be used to investigate noise, boundary effects and intrinsic alignments 
of galaxies (Catelan, Kamionkowski & Blandford, 2001; Crittenden, et al 2001; Croft & Metzler, 2001; 
Heavens, Refregier & Heymans, 2000), which will appear in both k and (3 modes. However, alignments 
appear to be a small effect at large redshift (Brown et al 2002, Heymans et al 2003). 

2 The Shear Power Spectrum 

2.1 Statistical properties 

We may define a shear covariance matrix by 

C ab (0) = (7«(O)7b(0)>, (8) 
where the indices (a, b) each take the values (1,2). Fourier transforming the shear field, 

7y (*) = jd 2 9 ll3 {0)e-* Le , (9) 

and decomposing it into k and j3 we may generate the shear power spectra: 

(k{£)k*{£')) = {2ir) 2 C^5 D {l-e.'), 

(n(£)f3*(£')) = (2ir) 2 Cf5 D (£-£'). (10) 



The parity invariance of weak lensing suggests that Cf' 3 = = 0. However other effects, such as 
noise and systematics, as well as intrinsic galaxy alignments may give rise to a non-zero Cf . The cross- 
correlation of k(£) and /3(£) is expected to be zero but it allows a second check on noise and systematics 
in the shear field. In particular finite field and boundary effects can lead to leakage of power between 
these three spectra. The shear power spectrum and the convergence power are related by Cj 1 — Cg K in 
the flat-sky approximation. For a spatially flat Universe, these arc in turn related to the matter power 
spectrum, P$(k,r) by the integral relation (see e.g. Bartelmann & Schneider 200f): 




where a is the expansion factor and r is comoving distance, rn is the comoving distance to the horizon: 
rjj = c J °° dz/H(z), where the Hubble parameter is given in terms of the matter density, Sl m , the vacuum 
energy density fly and the spatial curvature, CIk as H(z) = H a [(l + z) 3 Q m + (1 + z) 2 £Ik + ily] 1//2 . The 
weighting, W, is given in terms of the normalised source distribution, G(r)dr = p(z)dz: 

W{r)= { * dr' G(r') r —^. (12) 

J r ^ 

The covariances of the shear field, equation (8), are related to the power spectra by (Hu & White 2001) 

C u (r) = J-^L(crco S 2 2cpe + C^ sin 2 2 Vt -Cf sin 4^)| W{£)\ 2 e itf , 

C 22 (f) = J-0^(Crsm 2 2^ + C^cos 2 2^ e + Cfsm^ e )\W(e)\ 2 e d - f , 

Cl2(f) = /(0(^(C'r-C'f /3 )sin4^ + C'f cos4^)|TU(£)|V^, (13) 

where cos ipe = £-£x and £ x is a fiducial wavenumber projected along the x-axis. We have included here a 
smoothing or pixelisation window function; W(£) = jo(£ x dpix/2)jo(£ y 0pix/2), where jo = sin(x)/x is the 
zeroth order spherical Bessel function and 6> p i x is the smoothing/pixel scale. 

With these relations, the three shear power spectra C% K , Cf^ and can be estimated directly from 
shear data via a Maximum Likelihood approach (Hu & White 2001, Brown et al 2003). The dataset we 
have applied it to is the COMB017 survey. 

2.2 The COMB017 Dataset 

The data analysed is part of the COMBO-17 survey (Wolf et al. 2001), carried out with the Wide- 
Field Imager (WFI) at the MPG/ESO 2.2m telescope on La Silla, Chile. The survey currently consists 
of five 0.5° x 0.5° fields totaling 1.25 square degrees with observations taken in five broad-band filters 
(UBVRI) and 12 narrow-band filters ranging from 420 to 914 nm. The chosen filter set facilitates accurate 
photometric rcdshift estimation (a z w 0.05) reliable down to an i?-band magnitude of 24. During the 
observing runs the best seeing conditions were reserved for obtaining deep i?-band images of the five 
fields to allow accurate weak lensing studies. It is these i?-band images, along with the photometric 
redshift tables, that we make use of in this analysis. Gray et al (2002) discuss the procedure used to 
reduce the R band imaging data, which totalled 6.3 hours. The imcat shear analysis package was applied 
to our reduced image (see Gray et al 2002 for details). This resulted in a catalogue of galaxies with 
centroids and shear estimates throughout our field, corrected for the effects of PSF circularization and 
anisotropic smearing. We appended to this catalogue the photometric redshifts estimated for each galaxy 
from the standard COMBO-17 analysis of the full multicolour dataset. Wolf et al (2002) describe in 
detail the photometric redshift estimation methods used to obtain typical accuracies of a z = 0.05 for 
galaxies throughout < z < 1. 

2.3 The COMB017 shear power spectra 

To compare our data with other surveys we have estimated the shear correlation function, C(9), from 
COMB017 and compared these with the results of other surveys. The results are shown in Figure 1 
(LHS), where we see broad agreement with other surveys. 

Fig. 1 (RHS) shows the COMB017 shear power from a Maximum Likelihood analysis of the fields. 
Details of the analysis are given in Brown et al (2003). Three of our five band power measurements agree 
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Figure 1: LHS: The total correlation function, C(6), as measured from COMBO-17, along with the most 
recent cosmic shear measurements from the four other groups indicated. Beyond a scale of 1 arcmin, the 
measurements are in broad agreement. The correlation function predictions for a flat ACDM cosmology, 
for three values of the power spectrum normalisation (from top to bottom: as = 1-0, 0.9, 0.8) are also 
plotted. RHS: The cosmic shear power spectra from COMBO-17. Plotted on a linear-log scale are C^ K 
(circles), Cf (crosses) andCg (triangles) in 5 band-averaged band powers, as a function of multipole, £, 
estimated from the optimal combination of a maximum likelihood analysis of the four COMBO-17 fields, 
CDFS, SGP, FDF, and Sll. The error bars are estimated from the Fisher matrix. The solid curve is the 
shear power spectrum expected for a as — 0.8 normalised ACDM model. Once again, the C^ and C^ 
points have been slightly horizontally displaced for clarity (from Brown et al 2003) 



within their error bars with the as = 0.8 normalised ACDM model plotted. The high level of power 
measured at I ~ 4000 seems to be present in each of the fields we have analysed. The largest discrepancy 
between the model and the measurements occurs at t = 2000 where less than half the power is found. 
We can check our results for systematic effects by estimating the power in the [3-[3 modes, as well as the 
K-/3 cross correlation. In all of our band powers the (3-(3 correlation is below the detected signal in k 
and is consistent with zero in all but one band power at i ~ 4000. Similarly the n-f3 cross correlation 
is well below our measurement of shear power, and is consistent with zero except at I ~ 1000 where 
a significant detection appears. We conclude from the minimal power found in these spectra that our 
results are not strongly contaminated by systematic effects. Subsequently we have used the photometric 
redshift information to remove close galaxy pairs to remove any intrinsic alignment effects (Heyman et 
al 2003). This had a negligible effect on our results, indicating that for current shear surveys this is not 
a major source of systematic. 

2.3.1 Estimating TO and as 

We are now in a position to use our measured shear power spectrum estimates to obtain a joint mea- 
surement of the normalisation of the mass power spectrum erg, and the matter density f2 m , by fitting 
theoretical shear correlation functions and power spectra, calculated for particular values of these pa- 
rameters, to our measurements. Details are given in Brown et al (2003). Figure 2 shows the results 
of our parameter analysis in the as — fi m plane (dark thin solid lines) for COMB017. A best-fit gives 
a 8 = (0.72±0.09)(fi m /0.3)-°- 49 . 

2.3.2 Combination with the 2dFGRS and Pre-WMAP CMB experiments 

We can combine this confidence region with parameter estimations from other sources such as the 2dF 
Galaxy Redshift Survey (2dFGRS, Percival et al. 2002) and the various pre-WMAP CMB experiments 
(Lewis & Bridle 2002 & references therein). See Brown et al (2003) again for details of this combination. 
The results of this combination are shown in Figure 2 for an optical depth of r = 0.10. Note that 
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Figure 2: The likelihood surface of as and f2 m from combining the C'OMBO-17 dataset with the 2dFGRS 
and pre-WMAP CMB constraints. The dark thin solid contours are the constraints obtained from our 
shear power spectrum analysis described in the previous section. The lighter vertical contours are the 
constraints on VL m obtained from the 2dFGRS where we have applied the priors described in the text to 
the 2dF data. The dashed set of contours are the constraints from a compilation of six pre- WMAP CMB 
experiments where we have assumed an optical depth to reionization ofr — 0.10. The dark heavy contours 
are the 1,2 and 3a combined constraints from the three methods. The best fit values of Q m and as are 
also indicated (from Brown et al 2003). 

the 2dFGRS data we have used for this estimation constrains O m only and so the constraints shown 
on as come wholly from the cosmic shear and CMB measurements. We measure best-fit values of 
fi m = 0.27±o"i and a 8 = 0.73±°qI from the combined data, This compares well with the WMAP 
(+CBI, ACBAR, 2dFGRS and Lyman-a surveys) results of Q m = 0.27 ± 0.04 and a 8 = 0.84 ± 0.04. 
Although there is some overlap in these surveys there does seem to be good convergence in these results. 
However to push to the even high accuracies we can expect in future lensing surveys, the main sources 
of systematics must be identified and removed. 

3 Mapping the 3-D dark matter 

While the cosmic shear signal, with an accurate knowledge of the source redshift distribution, provides an 
accurate probe of the projected dark matter component and cosmological parameters, the combination 
of shear and source redshifts contains much more information. Wittman et al (2001, 2002) have demon- 
strated the utility of using source redshifts by inferring a cluster redshift from shear and photometric 
redshift information for galaxies in their sample. The importance of redshift information to remove in- 
trinsic galaxy alignments from shear studies has also been discussed by Heymans & Heavens (2002) and 
King & Schneider (2002a, b). Beyond this one would still like to image the full 3-D dark matter distri- 
bution over a large fraction of the Hubble volume. In particular this would aid seeing the cosmological 
growth of structure over cosmic time - something that can be calculated in the linear regime to high 
accuracy. In addition, while surveys such as WMAP have helped pin down cosmological parameter at 
high redshift, one would like to see them evolve. This is particularly important for parameters such as 
the equation of state of the universe, where its evolution can be used to distinguish between models for 
the dark energy. In fact it turns out that the 3-D mapping of dark matter is entirely feasible with good 
data. Here we outline the method proposed by Taylor (2001) and then an application to the COMB017 
data (see Taylor et al 2003). See also Hu & Keaton (2002) for a pixelised implementation to simulated 
data. 



Figure 3: Example of gravitational potential. Here we simulate the gravitational potential for two NFW 
clusters at r = (0.5,0.5,0.25) and (0.5,0.7,0.4); the upper panel is a 3-D representation of (y,z,$>) 
resulting from an x slice at x = 0.5°; the lower panel displays the $ values for (y,z) in greyscale. 

3.1 3-D reconstruction 

The inverse relation to equation (3) is 

$(r,r0) = ^d r r 2 d r (f>(r,r9), (14) 

which can be verified by substitution into equation (3) and integrating by parts. Here the Newtonian 
potential is evaluated at the source position, and I assume the lensing potential is smooth enough to 
allow differentiation by the radial derivative d r = 0.V. The shear field provides an estimate of the 
lensing potential up to an arbitrary quadratic function; 

0(r,r9) = J>(r,r9)+x(r,e), (15) 

where 

X (r, 9) = V(r) + n{r)0 x + p{r)9 y + v{r){6 2 x + 6 2 y ), (16) 

is a solution to 

(d i d j -±5Z#^ X (r,6) = 0, (17) 

and if), r], \i and v are arbitrary radial functions. This degeneracy is related to the so-called sheet-mass 
degeneracy which appears as a constant of integration, 4z/, when estimating the convergence from the 
shear field (Falco, Gorenstein & Shapiro 1985). However note that the degeneracy we find here will still 
appear, with v = 0, even if the convergence field is used to estimate the lensing potential. These terms 
can be removed by taking moments of the measured lens potential over the area of a survey. Defining 

Xp - q{r) = jJ A d20 ^ rd )( dP X + ° q y)> ( 18 ) 

where A is the area of a survey, an estimate of the lensing potential, with the mean, gradient and 
paraboloid contributions removed, is 

4> = $-x (19) 

where for a circular survey with radius 9, 

X = 4X0,0 (l - - ^ (l - |T) + J* (xi,oO x + XoM , (20) 

assuming that the true potential averages to zero. Hence 



$(r, rff) = \d r r 2 d r far, rO) - £(r, 0)) 



(21) 
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Figure 4: Top panel: reconstructed gravitational potential using the full shear field of Figure 6, in a 100 3 
grid. This is a (y, z) plane at x = 0.5°. Bottom panel: difference between input and recovered gravitational 
potential fields. 

is an estimate of the Newtonian gravitational potential, while 

S(r, r9) = 3^V 2 d r r 2 d r (£(r, r9)) ~ *(r, 9)) (22) 

is an estimate of the density field. The derivation of equations (21) and (22) are the main results of this 
paper, and demonstrate that the full, 3-D Newtonian potential and density fields can be reconstructed 
from weak lensing observations. In practice \ wu l n °t be a major problem for large surveys where the 
moments of the true lensing potential will average to zero. However for small fields this may not be so 
true, and one may instead wish to set the edge of the field to <f> = 0. 

3.2 Simulations 

Having set out the formal method for a 3-D reconstruction of the dark matter distribution, we now test 
it on simulations. Our simulation consists of a double cluster with a Navarro-Frenk- White profile (Figure 
3). Further details can be found in Bacon and Taylor (2003). 

3.3 Perfect dark matter reconstruction 

First we examine the accuracy of our method when the shear field is perfectly known everywhere. Figure 
4 displays the reconstructed $ and the difference between input and reconstructed fields. It can be seen 
that, with perfect knowledge of the shear field, a good reconstruction is achieved with our method. For 
the cluster at z = 0.25, the error is < 0.4% of the signal within a radius of 0.2 degrees of the cluster, 
except for the core pixel where the error is 11.5%. For the cluster at z — 0.5, the error is < 5% of the 
signal within a radius of 0.1 degrees of the cluster, except for the core pixel where the error is 17.6%. This 
is again due to the cusp at the cluster centre, which cannot be followed well by our averaged shear field. 
Nevertheless, it is clear that our method is successful in reconstructing cluster gravitational potentials in 
the absence of noise. 

3.4 Realistic Dark Matter Reconstructions 

Having demonstrated that the inversions of the shear field to obtain the gravitational potential is viable in 
the absence of noise, we now wish to add the two primary sources of noise present in lensing experiments: 
Poisson noise due to only sampling the field at a finite set of galaxy positions, and additional noise 
due to galaxies' non-zero ellipticities. As an example we use number densities for plausible space-based 
experiments (100 per sq arcmin), and an appropriate number of objects at each redshift slice. 

Figure 5 (LHS) shows the dark matter potential measurements for our notional space-based survey 
after Wiener filtering. We measure the z = 0.25 cluster with v = 4.2 at the peak (trough) pixel of its 



Figure 5: LHS Reconstructed dark matter potential for our space-based experiment, using galaxy prop- 
erties as in Figure 9, after Wiener filtering. Note the detections of clusters at (y, z) — (0.5, 0.2) and 
(0.7,0.4) with v = 4.2 and 2.1 respectively (crosses show trough minima). RHS Cross-section in redshift 
and angles of the 3-D dark matter potential for the COMB017 A901/2 super cluster field. Note that the 
lower large potential wells, at z — 0.2, correspond to the A901/2 supercluster, and the higher potential 
well, at z = 0.4, corresponds to an unknown background cluster. 

gravitational potential well, with v = 2.1 at the gravitational potential trough of the z = 0.4 cluster. 
We also find a substantial noise peak in the foreground at z = 0.1. While the recovery of the cluster 
gravitational potentials therefore constitutes a challenging measurement, we are indeed able to reconstruct 
useful information in the gravitational potential field itself. (The detection significance of these clusters 
is much higher than the measurement v at a given point in the cluster). 

3.5 Reconstructing the Dark Matter Potential in COMB017 

In this section, we will use these methods to calculate the dark matter potential for the COMB017 
supercluster A901/2 (See Gray et al 2002 for a 2-D shear analysis). For more details of the 3-D dark 
matter reconstruction see Taylor et al (2003). 

Fig. 5 (RHS) shows a cross-section through the final $ field. Note several significant features of this 
field. Firstly, we see that the gravitational potential trough associated with supercluster at z = 0.16 is 
clearly recovered at the bottom of the plot, with a peak pixel S/N of 2.9. Secondly, further analysis shows 
that the gravitational potential displays troughs at approximately the positions of the three component 
clusters of A901/2. Thirdly, we note that a second mass concentration is measured at z ~ 0.4, higher up 
in the plot. This has a peak pixel S/N of 2.8, so is certainly significant. 

While the signal-to- noise on this ground-based application of the 3-D reconstruction is low, this clearly 
demonstrates the method. With higher numbers of sources in space-based surveys, we can expect to map 
out the dark matter distribution in 3-D over a significant fraction of the Hubble volume. Such 3-D dark 
matter catalogues can then be used as the basis for constructing mass-selected cluster catalogues with 
no projection effects. 

4 Dark Energy from Lens Tomography 

While the combination of shear and redshift information can be use to directly map the 3-D dark matter, 
an orthogonal application can also be made to probe the dark energy component of the Universe. The 
dark energy makes itself felt via its effect on the evolution of the universe and in particular the Hubble 
parameter, H(a). This effects both the evolution of dark matter, and the geometry of the Universe. 
This effect has been utilised by, for example, Hu (2002), Huterer (2002), Abazajian & Dodelson (2002), 
Heavens (2003), Refregier et al (2003), Knox (2003) and Linder & Jenkins (2003). 



Recently Jain & Taylor (2003) have proposed a new method based purely on the geometric effect 
of dark energy on gravitational lensing. This involves a particularly simple cross-correlation statistic: 
the average tangential shear around massive foreground halos associated with galaxy groups and galaxy 
clusters. We show that cross-correlation tomography measures ratios of angular diameter distances over 
a range of redshifts. The distances are given by integrals of the expansion rate, which in turn depends 
on the equation of state of the dark energy. Thus the lensing tomography we propose can constrain the 
evolution of dark energy. Lens tomography has also been studied as a valuable means of introducing 
redshift information into shear power spectra as a general method for improving parameter estimation 
(see e.g. Seljak 1998, Hu 1999, 2002, Huterer 2002, King & Schneider 2002b). 

4.1 Method 

The dark energy has equation of state p — wp, with w = — 1 corresponding to a cosmological constant. 
The Hubble parameter H (a) is given by 



H(a) = H a 



n m a' 3 + n dc e- 3 Ii dlna ' {1+w{a,)) 



(23) 



where H is the Hubble parameter today. The comoving distance is r(a) = J^ a da' / a' 2 H (a') . 

We consider the lensing induced cross-correlation between massive foreground halos, which are traced 
by galaxies, and the tangential shear with respect to the halo center (denoted here 7): ujy.{8) = 
(8nf (0)7(0')) where nf(0) is the number density of foreground galaxies with mean redshift (zf), ob- 
served in the direction 9 in the sky and 5rif(0) = (rif(0) — fif)/hf. The angle between directions 9 and 
9' is 9. The cross-correlation is given by (Moessner & Jain 1998, Guzik & Seljak 2002): 



oj x (9) = 67r 2 a„ / H drW t {r)^^ f dk k P hm (r, k) J [kr9] , 
Jo a \ r ) Jo 



(24) 



where g{r) is the lensing geometry averaged over the normalized distribution of background galaxies 
W b ( X ) 

f™ r' -r 

g{r)=r dr'^^W h (r'), (25) 

J r ' 

Phm(r,k) is the halo-mass cross-power spectrum, and Wf is the foreground halo redshift distribution. 
The Bessel function J M has subscript fi = 2 for the tangential shear and /i = for the convergence. The 
measurement of the mean tangential shear around foreground galaxies is called galaxy-galaxy lensing. 
We will consider a generalization of this to massive halos that span galaxy groups and clusters. If the 
foreground sample has a narrow redshift distribution centered at r = rf, then we can take Wf to be 
a Dirac-delta function and evaluate the integral over r. All terms except g{rf,r^) are then functions 
of rf, the redshift of the lensing mass. The coupling of the foreground and background distributions 
is contained solely in .g(r f ,r b ). Hence if we take the ratio of the cross-correlation for two background 
populations with mean redshifts z\ and z 2 , we get 

^x,i(6>) = gijn) _ (n - r f )/n 

w x ,2(0) 92(ri) ~ (r 2 - n)/r 2 ' 

where wi,gi denote the values of the functions for the background population with mean redshift z\, 
and the second approximation is in the limit that the background galaxies also have a delta-function 
distribution. The above equations show that the change in the cross-correlation with background redshift 
does not depend on the galaxy-mass power spectrum, nor on 9. We can simply use measurements over a 
range of 9 to estimate the distance ratio of equation (26) for each pair of foreground-background redshifts. 
The distance ratio in turn depends on the cosmological parameters fide, w, and its evolution w' . 

4.2 Signal-to-noise estimate 

A simple way to estimate the signal-to-noise for the cross-correlation approach is to regard the foreground 
galaxies as providing a template for the shear fields of the background galaxies. For a perfect template (i.e. 
for high density of foreground galaxies and no biasing), the errors are solely due to the finite intrinsic 
cllipticities of background galaxies. Thus the fractional error in our measurement of the background 
shears is simply 

<fy <r e „„ „„_, /100\ 1/2 / 0.1 y /2 
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Figure 6: Contours in the fide — w plane for the fiducial tensing survey with / s k y = 0.1 . The inner contour 
assumes no evolution of dark energy, w a — 0, while the two outer contours marginalize over w a , with 
external constraints on Q v corresponding to cr(fl v ) — 0.01,0.03 (see text). The 68% confidence interval is 
shown in each of the contours. Contours in the w — w' plane for the fiducial lensing survey with / s k y = 0.1, 
as in Figure 6. The inner contour assumes Q v = 0.7, while the outer two contours marginalize over Q v 
as in Figure 6. Note that the parameter w a = 2 w' at z = 1 (from Jain & Taylor 2003). 



where the total number of background galaxies is A/total = n g A = n g / s ky^4 s k y , A is the survey area, 
the lensing induced rms shear (7) rms ~ 0.04, the intrinsic ellipticity dispersion o~ e = 0.3, and where the 
number density n g has units per square arcminute. Thus for the fiducial parameters / s ky = 0.1 and 
n g = 100, one can expect to measure the background shear to 0.1% accuracy at about 5-ct. We find that 
such a signal corresponds to changes in w of a few percent; hence this is the approximate sensitivity we 
expect in the absence of systematic errors. 



4.3 Dark energy parameters 

We perform a x 2 minimization over the mean shear amplitudes at the two background distributions for 
each foreground slice to fit for the dark energy parameters. The time dependence of w is parameterized 
as w — wq + w a (l — a) (Linder 2002). For comparison with other work we will compare w a to w' defined 
by w = w + w'z. Since the w' parameterization is unsuitable for the large redshift range we use, any 
comparison with it can be made only for a choice of redshift. For example, at z = 1, which is well probed 
by our method and is of interest in discriminating dark energy models (Linder 2002), w a = 2w' . For 
foreground slices labeled by index I and two background samples by 1 and 2, the \ 2 IS given by 



* 2 = E 



R°(zi,zi,z 2 y 2 



R(zi,zi,z 2 ) 



Uu (28) 



where R is the distance ratio of equation (26) for given values of Q v ,w and w a , and R° is the fiducial 
model with fl v — 0.7, w — —1, w a = 0. The weights Ui are given by 

1 2n 1 / i A( 7 )f 1 + 2n 2 .f i A( 7 )f 2 ' ^ 

where fi is the fraction of the survey area A covered by halo apertures in the i-th lens slice. The factor of 
2 in the denominator arises because we are using only one component of the measured ellipticity whereas 
of denotes the sum of the variances of both components. 

The results are shown in Figures 6. Figure 6 (LHS) shows the constraints in the Q v — w plane. The 
ellipses show the 68% confidence region given by A% 2 = 2.3. The elongated inner contour is for fixed 
w a = 0. The two outer contours marginalize over w a assuming external constraints on Q v from the CMB 
and other probes, corresponding to a(Q v ) — 0.01,0.03. 



Figure 6 (RHS) shows the constraints in the w — w a plane if fi„ is fixed, or marginalized over with 

— 1/2 

a(Cl v ) — 0.01,0.03. The corresponding accuracy on w and w can be scaled as f sky . For the case 

with a(Q v ) = 0.03, we obtain a(w) ~ 0.01/ s ^/ 2 and a(w a ) — 0.06/T.^/ 2 (at z = 1 this is equivalent to 

ct(w') ~ 0.03/ s ^/ 2 ). Note that the value of a on a parameter is given by projecting the A\ 2 — 1 contour 
on the parameter axis, which is smaller than the projection of the inner contour in Figure 6 (RHS) 
with A% 2 = 2.3. The scaling with / s k y in the parameter errors comes from the number of background 
galaxies (not sample variance), hence for given / s k y varying the depth of the survey scales the errors 

— 1/2 

roughly asn 9 . The results we have shown are for the fiducial rcdshift z = 0. A different choice of the 
fiducial rcdshift changes the relative accuracy on w and w a somewhat, because the degeneracy direction 
in the three parameters changes. A detailed exploration of different models of w(a) with finer bins in the 
background redshift distribution would be of interest. 

5 Summary 

In this paper I have presented three ways of combining weak lens shear data with source depth information. 
Firstly, one can use shear and redshifts to tie down the positions of source galaxies and remove this as 
a source of uncertainty in the measurement of the shear power spectrum. We have demonstrated this 
with an application to the COMB017 dataset (Brown et al 2003). Secondly, with source redshifts I have 
demonstrated on simulations and with the COMBO 17data, that the full 3-D dark matter distribution can 
be mapped (Taylor et al 2003). And finally with source redshifts the effect of the dark energy component 
of the universe can be explored via its effect on the geometry of gravitational lensing by taking ratios of 
cross-correlations of lens and source galaxies (Jain & Taylor 2003) . Clearly, from these applications and 
many more, the availability of source redshifts has opened up a whole new dimension for gravitational 
lensing. 
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